Semiclassical Foundation of Universality in Quantum Chaos 
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We sketch the semiclassical core of a proof of the so-called Bohigas-Giannoni-Schmit conjecture: 
A dynamical system with full classical chaos has a quantum energy spectrum with universal fluc- 
tuations on the scale of the mean level spacing. We show how in the semiclassical limit all system 
specific properties fade away, leaving only ergodicity, hyperbolicity, and combinatorics as agents de- 
termining the contributions of pairs of classical periodic orbits to the quantum spectral form factor. 
The small-time form factor is thus reproduced semiclassically. Bridges between classical orbits and 
(the non-linear sigma model of) quantum field theory are built by revealing the contributing orbit 
pairs as topologically equivalent to Feynman diagrams. 
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Fully chaotic dynamics enjoy ergodicity and thus visit 
everywhere in the accessible space with uniform likeli- 
hood, over long periods of time. Even long periodic orbits 
bring about such uniform coverage. Moreover, classical 
ergodicity provides quantum chaos with universal char- 
acteristics. 

Given chaos, quantum energy levels are correlated 
within local few-level clusters but become statistically in- 
dependent as their distance grows much larger than the 
mean level spacing A. The decay of correlations on the 
scale A is empirically found system independent, within 
universality classes distinguished by presence or absence 
of time-reversal (T) invariance 0, |2] ■ The corresponding 
universal long-time characteristics act on the Heisenberg 
scale Th = 27r?i/A, with h Planck's constant. 

Universal spectral fluctuations were conjectured as a 
manifestation of quantum chaos two decades ago 
Now, the semiclassical core of a proof can be given. Based 
on Gutzwiller's periodic-orbit theory 0, our progress 
comes with two surprises: one lies in its simplicity, the 
other in the appearance of interesting mathematics (non- 
trivial properties of permutations). Moreover, the often 
disputed intimate relation between periodic orbits and 
quantum field theory is confirmed for good. We thus 
expect the underlying ideas to radiate beyond spectral 
fluctuations, like to transport and localization. 

Technically speaking, we want to show that each com- 
pletely hyberbolic classical dynamics has a quantum en- 
ergy spectrum with the same fluctuations as a random- 
matrix caricature H^mt of its Hamiltonian, even though 
that caricature has nothing in common with the Hamilto- 
nian but symmetry (absence or presence of T invariance) . 
The theory of random matrices (RMT) J., devel- 
oped by Wigner and Dyson to account for fluctuations in 
nuclear spectra yields analytic results for correlators of 
the level density p{E), by averaging over suitable ensem- 
bles of random matrices. Simplest is the two-point corre- 
lator p{E)p{E') — p{E) p{E'), where the overlines denote 



ensemble average. Its Fourier transform with respect to 
the energy difference E — E', called spectral form factor 
K{t), is predicted by RMT for systems without time re- 
versal invariance (unitary class) and with that symmetry 
(orthogonal class) as 

ifu„i(T) = T , iforth(T) = 2t - T ln(l + 2t) , (1) 

respectively; here r is a time measured in units of the 
Heisenberg time Th ranging in < t < 1 . Note that the 
"orthogonal" form factor admits the expansion K{t) = 
2t — 2t^-|-2t^ . . . which converges for < r < i. Leaving 
larger times for future work we propose to show fidelity 
of individual chaotic dynamics to ^ . 

Of the many ways of doing RMT averages yielding 
the quantum field-theoretical non-linear sigma model 
deserves special mention since it yields a r expan- 
sion of the form factor equivalent to the semiclassical 
expansion to be developed here. The model points to 
analogies between hyperbolic dynamics and the motion 
of electrons in disordered media. In fact, the equivalence 
of semiclassical and field-theoretic expansions was first 
suggested in the context of disordered metals (gj]- 

Our starting point is Gutzwiller's representation of the 
level density of a hyperbolic system by a sum over its 
classical periodic orbits, p{E) oc Re^^ Aj&^^^^ , with 
Sy{E) the action and the (dimensionless) stability 
amplitude of the 7th orbit. The form factor K{t) is the 
double sum 
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with Ty{E) the period of 7; the angular brackets demand 
averages over the energy and over a time interval AT <C 
Th- We aim at evaluating the periodic-orbit sum |(2Jl in 
the semiclassical limit ?i — > 0, — *■ 00 with Tj/Th — 
const. That limit and the averages indicated eliminate 
noise due to orbits with I I ^ h and purge K{t) 

of system specific features. 
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For the formal double sum in ^ to converge to the 
RMT prediction it must be structured into contri- 
butions from families of orbit pairs, such that each term 
of the T expansion of K{t) comes from a specific set of 
families. The simplest family contains the diagonal pairs 
{7, 7} and, given time reversal invariance, {7, Tj}, where 
T7 is the time reverse of 7; it yields Berry's (9( "diago- 
nal approximation" K^]} — kt where k = 1 without and 
K = 2 with T invariance, due to the doubling of pairs in 
the latter case. It is here, when summing over the "diag- 
onal pairs" that we first meet the ergodicity of long peri- 
odic orbits, through Hannay's and Ozorio de Almeida's 
(HOdA) [13 sum rule {J^^ \A^\^6{t - ^)) = t. In a 
paradigmatic breakthrough, Sieber and Richter gave 
the family responsible for the term of T invariant dy- 
namics; it is on the basis of their insight that we could 
find and account for all other families. 

We first turn to the unitary class and propose to 
demonstrate that all families of orbit pairs individually 
contributing to higher orders r" collectively cancel for 
n > 1. To ease our task we assume two freedoms. 

Long orbits have lots of close self-encounters. We speak 
of an Z-encounter when I orbit stretches get and stay close 
for as long as their chaotic, i.e. exponential divergence 
permits (see Fig. 1). Since the closest approaches dis- 
cernible quantum mechanically have an action scale h we 
expect relevant encounter durations tcnc of the order of 
the Ehrenfest time Te In i, with A the Lyapounov 

rate of divergence. Departing from and ending on the 
21 "ports" of an Z-encounter are / "loops" with durations 
of the order of the period T and thus of the Heisenberg 
time Th = 2nfi/ A. = Q,/2TTh, where is the volume of 
the energy shell. Different encounters must be considered 
as separate: overlap of any two would yield a single one 
with more internal stretches. More generally, an orbit 
must leave an encounter before reentering it or another 
one after traversing an outside loop. 

Self-encounters lead us from an orbit 7 to partners 7'. 
Two orbits in a pair {7, 7'} are practically indistinguish- 
able in the loops outside encounters; they only differ 
within comparatively short encounters, by their connec- 
tions of the outside loops. The action difference — S^i 
can thus be of order Ti. ReshufHing intraencounter con- 
nections of 7 either yields a partner orbit 7' or a pseudo- 




FIG. 1: A triple encounter = 3) in the energy shell and 
its Poincare section for an orbit pair 7,7'. Inset: Global 
appearance of the pair and the generating encounter. 



orbit decomposing into shorter orbits; pseudo-orbits. 
Fig. 2, are not admitted to the Gutzwiller sum (0). 

Calling vi the number of Z-encounters {I > 2) within 
which 7 and 7' differ by connections of the coinciding 
outside loops we write V = vi for the total number 
of such encounters and L = Ivi for the number of or- 
bit stretches within encounters (equalling the number of 
loops outside). We shall see that the families of orbit 
pairs with fixed n = L — V + 1 exclusively contribute to 
t" in K{t). To calculate those contributions and check 
that they sum up to zero for n > 1 in the absence of 
T invariance, first we must understand the phase-space 
structure of self-encounters and, second, the combina- 
torics of counting proper partners must be mastered. 

We begin with a closer analysis of self-encounters ■ 
Drawing a Poincare section, two dimensional for two- 
freedom systems, through an Z-encounter we see the I 
orbit stretches of, say, 7 pierce through that section in 
I points Xi = {ui,Si), i = 1...Z; one of these can be 
chosen as the origin of coordinate axes spanned by the 
unstable and stable manifolds of 7 through Xi = (0,0). 
For an encounter to be close we require |si| < c, with 
the bound c small enough for the motion along the I orbit 
stretches to allow for mutually linearized treatment. The 
I — 1 piercings Xi ^ Q uniquely determine (i) the duration 
tclc(u, s) of an Z-encounter as a (logarithmic) function of 
the Ui, Si [13, (n) the piercings of the partner orbit (s), 
and (iii) the contribution to the action difference 
There is a canonical transformation u, s H, s diagonal- 
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izing the action difference to AS", 

tcnc and AS'onc are canonical invariants. 

We may characterize a set of encounters by a vector 
1/ whose components are the numbers vi of /-encounters. 
We define a (weighted) number w{u, s)d'"~^ud^~^s of 
encounter sets with fixed vector v and temporal order of 
the L visits of the V encounters inside an orbit of period 
T; it contains a factor for each encounter involved, the 
fraction of its duration which the corresponding unstable 
and stable components spend in the intervals [ui , Ui + 
dui] and [si,Si -I- dsi]. That number is determined by 
ergodicity as follows. A piercing of 7 through a section 
will be found with the uniform probability dtiduds/Cl in 
a time interval [ti , ti + dti] and in the area element [u,u + 
du] X [s, s + ds]. We integrate the product oi L — V such 
probabilities over the L times ti, here ti G [0,T] while 
we restrict the L — 1 other ti to (i) a specific order in the 
interval [ti, ti+T] and (ii) by minimal separations (due to 
the ban of encounter overlap). To get the dimensionless 




FIG. 2: Illustration of the counting problem: pseudo-orbits 
(here, dashed) must be eradicated. 
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weight w{u, s)d^ ^ud'" we divide the L-fold integral 
by the product J| tone of durations of the V encounters 



T 



w{u, s) 



L-l 



one 



(3) 



here, the restrictions mentioned in effect reduce the inter- 
val accessible to L — 1 integration variables by the cumu- 
lative duration toxci — '^onc of the L intra-encounter 
stretches. We note that the contribution Itcnc of each l- 
encounter to texci depends on the I ~ 1 points of piercing 
X = (m, s) 7^ of 7 through the pertinent section. We 
shall see below that the non-vanishing duration tcxci is, 
even though a small correction to the period T in s), 
of decisive importance for spectral universality. 

Remarkably, s) results as independent of the order 
of the L passages of 7 through the V self-encounters. We 
can therefore proceed to the number of self-encounters 
with fixed v and u, s irrespective of the order of visits, 
N{v)w{u^ s), by accounting for a multiplicity N{v). 

The number N{v) brings up a combinatorial problem 
with a shade of topology mixed in (partner orbits must 
be connected), decoupled from the phase-space consider- 
ations yielding w{u, s). When stating that 7 and 7' differ 
in vi /-encounters, I — 2,3, . . . we leave open (i) the order 
in which the L encounter stretches are passed (in partic- 
ular, which is the first) and (ii) how the intraencountcr 
connections of 7 are reshuffled in 7'. The number of 
possibilities left is N{v); we can determine it by running 
through all different orderings of visits as well as through 
all intraencounter connections other than the one realized 
by 7 and checking, in each case, whether an orbit or a 
pseudo-orbit results. For only a few vi non-zero, N{v) is 
easily found with paper, pencil and patience. For gen- 
eral V, the permutation problem at issue can be attacked 
recursively as discussed in the technical note below. 

With the help of the weight N{v)w{u, s) we re- 
place the sum over orbit pairs in ^ as J^-yj' 
T,-yT,v^i'^)Id^~^ud^~^sw{u,s)L-\ The number of 
loops L had to be divided out here, since the L choices of 
one intra-encounter stretch as the first yield the same 
partner orbit. The summand simplifies as A^A*, — > 
A-yl^, T.y+T.y' ^ 2Tj siucB in contrast to the action dif- 
ference AS* the prefactors and periods suffer no relative 
discrimination by a small quantum unit. We may also 
invoke the HOdA sum rule already met above, to do the 
sum over 7, and thus find the form factor as 



K{t)-t 



^Y.N{v)rd^-''ud'^-^i 
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Only a single term of the multinomial expansion of 
(T — texci)^"^ in w{u, s) survives the limit h ^ 0, the one 
which cancels the denominator Y[ tine and comes with the 



factor (T/r2)^~^; all other terms vanish, either because 
they involve too low orders in the period T and thus ex- 
tra factors h besides a power of T/Th or because they 
oscillate rapidly and are annulled by averaging over a 
small time window. The remaining integral becomes ele- 
mentary after the canonical transformation diagonalizing 
A^, f_J^-^ud^-^se'^S/^ = (?i/'^dudse^"^/'"')^"^^ 
{2TTh)^'^ . Writing {T2nh/n)^-^ = t^"^ we get the 
series K{t) — t + X]^2 ^"T" with the coefficient 
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{n=L-V+l) 



(n-2)! 



governed by ergodicity, combinatorics and topology, 
given the existence of sets of separated close self- 
encounters. 

The vanishing of the foregoing sum over families of 
partner orbits is a property of the permutation group 
which to the best of our knowledge was never noticed 
before. We sketch the surprisingly simple proof of Kn = 
for n > 0, based on a recursion scheme for N{v), in the 
technical note below. Universal spectral fluctuations are 
thus established for dynamics without T invariance. 

The orthogonal class of T invariant dynamics can be 
treated similarly. We must generalize the notion of self- 
encounters to include orbit stretches close to the others 
up to time reversal. Configuration-space pictures prove 
useful: an orbit stretch may be depicted by an arrow 
— ' . While all /-encounters admitted in the unitary case 
involve / parallel such arrows (like =£: or 3) we now 
face, in addition, arrows with opposite directions (like 
:;=i: or ^). Likewise, loops in between self-encounters 
of an orbit 7 appear nearly unchanged in partner (s) 7', 
except that the senses of traversal may be opposite. 

The multiplicity N{v) of encounter "classes" v leads 
to a permutation problem slightly more complicated 
than in the unitary case. Again, all classes with fixed 
n = L - V + 1 contribute to t" . The results (gj and © 
reappear with an additional factor 2 due to the fact that 
with 7' a partner so is T7'. As discussed in the technical 
note, a recursion relation arises for Kn which yields the 
random-matrix form factor for the orthogonal universal- 
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We would like to 



ity class, iCorth = 2r-|-X;„>2 „-i 
underscore that in establishing both the unitary and the 
orthogonal form factor as universal we have accounted 
for all orbit pairs whose members differ by nothing but 
the way almost coinciding (up to time reversal) loops are 
connected within close self-encounters. 

The r expansion of iCorth converges for < r < i. The 
summed up logarithm remains valid, by analytic contin- 
uation, up to the next singularity. Neither the locus of 
that singularity (r = 1) nor the form factor for r > 1 
can be found within the r expansion. We underscore 
once more ergodicity and hyperbolicity as our basic as- 
sumptions; in addition, strong action degeneracies as for 
dynamics with Hecke symmetries must be excluded 13]. 
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We must discuss the relation of our semiclassical work 
to the so-cahed zero dimensional sigma model of quantum 
field theory 0,H,I3- The relevance of the sigma model for 
us lies in similarities of its perturbative implementation 
to our semiclassical expansion. A perturbative evaluation 
of the sigma model involves Wick's theorem which can 
be shown to entail a recursive reduction scheme equiva- 
lent to the topological and combinatorial problem yield- 
ing our multiplicity N{v). Moreover, our orbit pairs cor- 
respond to the Feynman diagrams depicting terms of the 
perturbative treatment of the sigma model, with our l- 
encounters and the outside loops the analogs of vertices 
(with 21 ports) and propagator lines, respectively. Order 
by order in r, our families of orbit pairs are equivalent to 
the Feynman diagrams of the sigma model 

Technical note: We want to set up the permutation 
problem yielding the mutiplicity N(v), first for the uni- 
tary case. To that end, starting from an arbitrary orbit 
stretch in some encounter we number the L stretches in 
the order of visits by 7. More precisely, we denote en- 
trance ports of encounters by 1,2, ... ,L and exit ports 
by 1', 2' , . . . , L' , such that the k-ih stretch of 7 connects 
ports k and k' . In a partner orbit 7' of 7 port k is con- 
nected to a port ^ k' . The L intra-encounter stretches 
of 7 thus correspond to the trivial permutation P'^ = 

(i' 2' L') '^hile a partner 7' will have intra-encounter 
connections according to Pe — ) Pe- Since re- 

connections take place only within encounters, Pe must 
be composed of vi cycles of length I, I ~ 2,3 . . .. (Math- 
ematically speaking, Pe must belong to the conjugacy 
class 2^23"^ . . .1""' of the group of permutations of L ob- 
jects corresponding to the cycles defined by the vector 
v.) 

The loops common to 7 and its partners 7' are associ- 
ated with the permutation Pl ^ \ \ " \ ^ ■ The whole 

of 7 is represented by the product PePl = (231) 
and that product is a single-cycle permutation since 7 is 
a single orbit, rather than a decomposing pseudo-orbit. 
Moreover, the product PePl describes a connected part- 
ner 7' rather than a decomposing pseudo-orbit if and 
only if it is single-cycle as well. The multiplicity N{v) is 
thus found by running Pe through all possibilities and 
counting only those for which PePl is single-cycle. 

In the orthogonal case there are 2'~^ distinct orienta- 
tions of the stretches of an ^-encounter. After combining 
oriented loops with reconnections of encounter stretches 
we must again check connectivity. 

We have established a general recursion relation for the 
multiplicity N{v), both in the unitary and the orthogonal 
case, by following the change of N{v) as (i) two encoun- 
ters unite {vi vi-l,vi> —>■ vi> - l,vi+i'-i + 
in short v — > t/^''' 1), (ii) as an encounter splits into two 
and (iii) as an Z-encounter becomes an {I — l)-encounter 
by uniting two of its orbit stretches. The relations are 



best written for N{v) = N{v){-1)^ 1]; [^(" ~ 2)!]-i; 
note that = -^=""1) N{v). 

In the unitary case we only need the special vari- 
ant of the general recursion relation concerning 2- 
encounters merging with Z-encounters, where the recur- 
sion reads V2N{v) + J2i>2 K^i+i + l)N{v^^^'^^) = 0. Sum- 
ming over V with fixed n = L — V + I we obtain 
^("=^-^+1) [v^Niv) + Y.i>A^i+i + i)^(t/[''^i)] - 0. 
In the foregoing double sum we can replace 
and thus vi+i -|- 1 — > The sum now runs over all 

V with > 0; however, the latter restriction is imma- 
terial due to the factor in the summand. We thus 

obtain [^■2 + Y.i>2 lvi+i\N{v) = 0; here, the 

term in the square bracket equals n — 1. The resulting 
identity (n — \)Kn — implies Kn = for n > 1. 

In the orthogonal case we need two special cases of the 
general recursion relation, to account separately for "dis- 
appearance" of 2-encounters and 3-encounters. A suit- 
able linear combination yields Kn+i = — ^^"^^^^ K^, the 
latter recursion gives the random-matrix form factor. 
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